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The relative entropy of entanglement is defined in terms of the relative entropy between an 
entangled state and its closest separable state (CSS). Given a multipartite-state on the boundary of 
the set of separable states, we find a closed formula for all the entangled states for which this state 
is a CSS. Our formula holds for multipartite states in all dimensions. For the bipartite case of two 
qubits our formula reduce to the one given in Phys. Rev. A 78, 032310 (2008). 
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I. INTRODUCTION 

Immediately with the emergence of quantum informa- 
tion science (QIS), entanglement was recognized as the 
key resource for many tasks such as teleportation, su- 
per dense coding and more recently measurement based 
quantum computation (for review, see e.g. [1,,2^). This 
recognition sparked an enormous stream of work in an ef- 
fort to quantify entanglement in both bipartite and mul- 
tipartite settings. Despite the huge effort, except the 
negativity [31j (and the logarithmic negativity [4J) closed 
formulas for the calculation of different measures of en- 
tanglement exist only in two qubits systems and, to our 
knowledge, only for the entanglement of formation [5]. 
Moreover, the discovery that several measures of entan- 
glement and quantum channel capacities are not addi- 
tive [6l [7] , made it clear that formulas in lower dimen- 
sional systems, in general, can not be used to determine 
the asymptotic rates of different quantum information 
tasks. Hence, formulas in higher dimensional systems 
are quite essential for the development of QIS. 

Among the different measures of entanglement, the rel- 
ative entropy of entanglement (REE) is of a particular 
importance. The REE is defined by [8,: 



En{p) = mmS{p\\a') 



S{p\\a) 



(1) 



where V is the set of separable states or posi- 
tive partial transpose (PPT) states, and S{p\\cf) = 
Tr (plogp — p log cr). It quantifies to what extent a given 
state can be operationally distinguished from the clos- 
est state which is either separable or has a positive par- 
tial transpose (PPT). Besides of being an entanglement 
monotone it also has nice properties such as being asymp- 
totically continuous. The importance of the REE comes 
from the fact that its asymptotic version provides the 



unique rate for reversible transformations [9 . This prop- 
erty was demonstrated recently with the discovery that 
the regularized REE is the unique function that quantify 
the rate of interconversion between states in a reversible 
theory of entanglement, where all types of non-entangling 
operations are allowed [10 . 

The state a = cf{p) in Eq. (fTl) is called the closest sep- 
arable state (CSS) or the closest PPT state. Recently, 
the inverse problem to the long standing problem pTj of 
finding the formula for the CSS a{p) was solved in [12] 
for the case of two qubits. In [12 the authors found 
a closed formula for the inverse problem. That is, for 
a given state < a on the boundary of 2-qubits sepa- 
rable states, dV, the authors found an explicit formula 
describing all entangled states for which a is the CSS. 
Quite astonishingly, we show here that this inverse prob- 
lem can be solved analytically not only for the case of two 
qubits, but in fact in all dimensions and for any number 
of parties. 

We now describe briefly this formula. Denote by H^ 
the Hilbert space of n x n hermitian matrices, where the 
inner product of X, F G H^ is given by Tr XY. Denote by 
Hn,+,i C iln,-\- C Hn the convex set of positive hermitian 
matrices of trace one, and the cone of positive hermitian 
matrices, respectively. Here n = nin2 • • • n^ so that the 
multi-partite density matrix p G H^^+^i can be viewed as 
acting on the s-parties Hilbert space C^^ ^C^^ (8)- • •(8)C^^ . 

Let < cr G H^^+ (i.e. a is full rank). Then for any 
<j' G H^ and a small real e we have the Taylor expansion 
of 

log(cr + ea) = logcr + eL^{a) + 0{e'^). 

Here L^r : H^ ^ H^ is a self- adjoint operator (defined 
in the next section), which is invertible, and satisfies 

Assume now that < a G dV (later we will extend 
the results for all a G dV; i.e. not necessarily full rank). 
Then, from the supporting hyperplane theorem, a has at 
least one supporting hyperplane, ^ G H^, of the following 
form: 
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Tr(0(j') > Tr(0(j) = V a' eV 



(2) 



where (j) is normalized; i.e. Tr^ = 1. For each such ^, 
we define the family of all entangled states, p{x^a)^ for 
which (J is the CSS: 



p{x^a) = a — xL^^{(j))^ < x < x^ 



(3) 



Here, Xmax is defined such that p{x^s,^^a) G H^^+^i and 
p(^max, cr) has at least one zero eigenvalue. We also note 
that TtL~^{(J)) = 0. Moreover, for the case of two qubits, 
(p is unique and is given by (j) = {\(p){(p\) , where T is 
the partial transpose, and \(p) is the unique normalized 
state that satisfies (J^\(p) = 0. Hence, for the case of two 
qubits our formula is reduced to the one given in [12] , by 
recognizing that for this case, the self-adjoint operator 
L~^ is given by the function G{a) of Ref. [12 . 

This paper is organized as follows. In the next section 
we discuss the definition of L^r. In section HI we find 
necessary and sufficient conditions for the CSS and in 
section IV we prove the main result for the case were the 
CSS is full rank. To illustrate how the formula can be 
applied, in section V we discuss the qubit-qudit 2 x m 
case. In section VI we discuss how to apply the formula 
for tensor products. In section VII we discuss the singular 
case, and show that the CSS state can be described in a 
similar way to the non-singular case. We end in section 
VIII with conclusions. 



II. DEFINITION OF L^ 

Let < a G Hn,+ . Fix /3 G H^. Let t G {-e,e) for 
some small e = e{a) > 0. Rellich's theorem, , e.g. [15], 
yields that log(a + 1/3) is analytic for t G (—£,£). So 



log(a + t(3) = log a + tLa{(3) + O(t^). 



(4) 



Here Lq, : H^ -^ H^ is the following linear operator. In 
the eigenbasis of a, a = diag(ai, . . . , a^) is a diagonal 
matrix, where ai, . . . , a^ > 0. Then for /3 = [bij]fj=i we 
have that 

[^a[P)\ki = hi , /c,/ = l,...,n. (5) 

ak — ai 

Here we assume that for a positive a, ^^^^Z^^^^ = 
a. 



a ' log a— log a 

Equivalently, for a real diagonal matrix a = 
(ai, . . . , a^) > define the real symmetric matrix 



[^(^)]m=i = 



log ak - log ai 

ak - ai 

ak - CLi 
log ak - log ai ' 



(6) 



Then, La{/3) = /3 o T{a), where /3 o r] is the entrywise 
product of two matrices, sometimes called the Hadamard 
product of matrices. Note that L^ is an invertible oper- 
ator, where L~^{l3) = /3 o S{a). 




FIG. 1: A schematic diagram of separable states (blue) and 
entangled states (orange). Most points on the boundary, like 
the points D and F, have a unique supporting hyperplane 
(which is also the tangent plane). The point F is the CSS of 
all the points on the red line. Some of the points, like the 
points C and E, have more than one supporting hyperplane. 
The point E is the CCS of all the points in the shaded green 
area. Some points on the boundary, like the points A and B, 
can not be a CSS; for example, separable states of rank 1 (i.e. 
product states) are on the boundary of separable states, but 
can never be the CSS of some entangled state. 



III. A NECESSARY AND SUFFICIENT 
CONDITION FOR a{p) 

We start with a necessary and sufficient condition the 
CSS, cr(p), must satisfy. 

Theorem 1. Let < p e Hn,+,i\P. The state < 
a{p) eV is a solution to Eq. ^, if and only if a = a{p) 
satisfies 



maxTrcr'Lcr(p) 



TraL,(p) = l. 



(7) 



Remark. We will see later that the assumptions that < 
p and < (T are not necessary. 

Proof First, note that Lcr{o-) = J^, and L^j is a self- 
adjoint operator. Hence, TrcrLcr(p) = TrLcr(cr)p = 
Tr(p) = 1. Now, let a^ G V. Since P is a convex set, it fol- 
lows that for every t G [0, 1], {l-t)a^ta' = a^t{a'-a) G 
V. Thus, applying Rellich's theorem for a small t > 
gives 

log(a + t{a' - a)) = log a + tL,{a' - a) + 0{t^). (8) 

If cr is a solution to Eq.(IT]), we must have Trploga > 
Trplog [cr + t(cr' — cr)] which together with Eq.Q im- 
plies that for a small positive t, tTi pL^{g' — cr) < 0. 
Dividing by t gives Ti pL(j{o-') < TrpLcr(cr) = 1. This 
completes the necessary part of the proof since L^r is self- 
adjoint. The sufficient part of the proof follows directly 
from the construction of p(x, a) in Theorem Is] D 

The proposition above leads to the following intuitive 
corollary: 

Corollary 2. Let < a be a CSS of an entangled state 
p. Then, a G dV. 



Proof. Assume that the CSS a is an interior point of V. 
Thus, for each a' separable (1 — t)a-\-ta^ is separable for 
all small |t|, where t is either positive or negative. Hence, 
instead of Eq.Q we get the identity Tt L(j{p)a' = 1 for 
all separable states a' G V. This yields that Lcr(p) = In- 
Hence p = a which is impossible since p was assumed not 
to be separable. D 



IV. MAIN THEOREM 

In the following we prove the main theorem for the 
case where the CSS, a^ is full rank. Note that if p is full 
rank then a also must be full rank. 

Theorem 3. (a) Let < a G dV, and let p G H^^+^i. 

Then, 

Er{p) = S{p\\a) {i.e. a is the CSS of p) 
if and only if p = p{x,a), where p{x^a) is defined in 

(b) If p > (i.e. full rank) than the CSS is unique. 

Proof, (a) We first assume that a is a CSS of p. Re- 
call that any linear functional ^ on H^ is of the form 
^{X) = Tr(0X) for some (j) e H^. Now, since V is 
a closed convex subset of Hn,+,i it follows (from the 
supporting hyperplane theorem) that for each boundary 
point a G dV there exists a nonzero linear functional 
on <l> : H^ ^ R, represented by ^ G H^, satisfying the 
following condition: 



<l>(cr) < ^(a ) for ah a' G V. 



(9) 



Note that the equation above holds true if (j) is replaced 
by — a/ (this is because Tra = Tra' = 1). Moreover, 
since 7^ we can normalize it. Therefore, there exists 
(/) G H^ satisfying ^ and the normalization 



Tr( 



1. 



(10) 



For most a on the boundary dV, the supporting hyper- 
plane of P at cr is unique (see Fig. 1). This is equivalent 
to say that ^ G H^ satisfying the conditions in Eq. ^ 
and (10) is unique. However, for some special boundary 
points a G dV^ there is a cone of such (j) of dimension 
greater than one satisfying Q (see Fig. 1). 

Now, denote (j)^ := —{La\p) — I). Since we assume 
that a is a CSS of p we get from Eq. (JTl the condition 
Tr((/)VO > Tr((/)V) = 0. Recah that LJ^cr) = I. Hence 
(j)^ = —L(j{p — a). Since p ^ a and L^r is invertible, it 
follows that 0' 7^ 0. Hence, (j)' can be normalized such 
that 0' = x^, where (j) satisfies Eq. ( 10 ) and x > 0. Apply 
L-i to ^ to deduce ^. We remark that TrL-i(^) = 0. 
Indeed 

= Tv4>a = Tv4>L-\l) = TtL-\<P)I = TV i^^ (</.). 

Assume now that < a G dV^ and let ^ be a sup- 
porting hyperplane at a, satisfying Eq. Q and (10). Set 



p = p{x^o-) as in Eq.([3|. We want to show that for this 
p, Er{p) = S{p\\a). Recall first that the relative entropy 
S{p\\o-') := Tr (p log p) — Tr(p log cr') is jointly convex in its 
arguments [16, Thm 11.12]. By fixing the first variable 
p we deduce that Tr(plogcr') is concave on V. Consider 
the function 

fit) := Tr(plog((l - t)a + ta')), t G [0, 1], 

where p > is given by Eq. (|3| and x > 0. The joint con- 
vexity of the relative entropy implies that f{t) is concave. 
Now, to show that the minimum of S{p\\o-') is obtained 
at cr' = cr, it is enough to show that /(O) > /(I) for each 
a^ G V. To see that, we first show show that /'(O) < 0, 
which then, combined with concavity of /, implies that 
/(O) > /(I). For small t we have 

log(cr + t{a - a')) = log cr + tL^{a' - a) + 0{t^). 

Hence, 

/(0) = Tr(pL.(a'-a)) = 

= Tr [L,{p){a^ - a)] = Tr [(/ - x^){a^ - a)] 
= X Tr(^cr - (j)a') = -x Tr(0crO < 0. (11) 

This completes the proof of part (a). Moreover, Er{p) = 
Tr(plogp) — Tr(p log cr) > 0, since p ^ a. Hence p is 
entangled. 

(b) This part follows from the strong concavity of log a 
(see appendix [A]) . Therefore, from Corollary 11 of ap- 



pendix A, it follows that for a fixed entangled state p > 0, 
the function Tr p log a is a strict concave function on the 
open set of all strictly positive Hermitian matrices in H^. 
Hence, if both a and a^ are CSS of p, then both are full 
rank and we have Trploga = Tiploga'. Hence, for 
t G (0, 1) we set a" = ta -\- {1 — t)a' and from the strong 
concavity 

Trploga' > tTrplogcr + (l-t)Trploga =Trplogcr , 

in contradiction with the assumption that a is a CSS. D 

Corollary 4. Let < p G Hn,+,i be entangled state 
and let a be a CSS of p. Then, a is also the CSS of 
p{t) = tp + (1 - t)cr for all t G [0, tmax]; where tmax > 1 
is the maximum t such that p{t) > 0. 

Proof. Since a is the CSS of p, from theorem [3] we have 
p = p(x, a) for some x. Hence p{t) = p(tx, a) is of the 
same form. From theorem [s] a is a CSS of p{t). D 

Remark. A weaker version of the corollary above was 
proved in [8]; note that here t can be greater than one. 



V. BIPARTITE PARTIAL TRANSPOSE 

ANALYSIS 

In the following we show how to apply Theorem [S] to 
specific examples. In particular, we focus on the bipartite 



case (i.e. n = 711122) and we will assume that V in Eq.Q 
is the set of PPT states. In the 2x2 and 2x3 case, 
V is also the set of separable states [17]. The boundary 
of the PPT states is simple to characterize. If a G V 
satisfies a > and also a^ > 0, where F is the partial 
transpose, then a must be an interior point of V. If on 
the other hand a or a"^ are singular, then a must be on 
the boundary of V. We therefore have: 

dV= (a eV det{a^a) = o| . 

Suppose now that < a G dV. Hence, a^ has at least 
one zero eigenvalue. Let |(^) be a normalized eigenstate 
corresponding to an eigenvalue zero and define an Her- 
mitian matrix (j) = {\(p) {(p\)'^ . Since the partial trans- 
pose is self- adjoint with respect to the inner product 
(/), p') = Tr(pp'), it follows that (j) satisfies Eq.([2]) and 
is normalized (i.e. Ti (j)'^ = 1). That is, (j) represents the 
supporting hyperplane at a. Note that if a^ has more 
than one zero eigenvalue than clearly (p is not unique and 
in fact there is a cone of supporting hyperplanes of V at 
a (see points C and E in Fig.l). To illustrate this point 
in more details, we discuss now the case where ni = 2 
(i.e. the first system is a qubit) and n2 = m. 

In the 2 X m case, we can write any state a G H2m,+,i 
using the block representation of 



a ■■ 



A B 
B^ C 



^(^(2m)x(2m)^ A, 5, C G 



(12) 



and A^ = A, C^ = C. The partial transpose of a is 
given by (here the partial transpose corresponds to the 
transpose on the first qubit system; i.e. it is the left par- 

A B^ 

B C 



tial transpose): <j^ 



. The following theorem 



shows that cr^ can have more then one zero eigenvalue. 

Theorem 5. Let m > 2. If a > and a^ > 
then rank a^ > m -\- 1. Furthermore^ for each k = 
0, . . . , m — 1 there exist strictly positive hermitian ma- 
trices a G Il2m,+,i such that a^ > 0, rank a^ = 2m — k. 

Proof. Recall that since a G H2m,+,i is strictly positive 
definite we have A > 0. Hence, <j and <j^ are equivalent 
to the following block diagonal hermitian matrices 



A 
C-B^A- 


-^B 






I 0' 


' A B' 
B^ C 




A 

C-BA-^B\ 




I 0' 
-BA-^ I 


A 
B 


5t- 
C 





I 

-B^A- 




1 / 



it 



respectively. Hence 



cr>0 < 
cr^ >0 



/ 
-BA-^ I 



C -B^A-^B>{) 
> C -BA-^B^ >0 



Note first that C ^ BA~^B^ . Otherwise, we get that 
BA-^B^ > B^A-^B, and since B^A'^B > it follows 
that each eigenvalue of BA~^B'^ must be positive and 
the i — th eigenvalue of BA~^B^ must be strictly greater 
then the i — th eigenvalue of B^A~^B. This can not be 
true since det BA~^B^ = det B^A~^B. Hence rank a^ > 
771 + 1. This complete the first part of the theorem. 

Next, let ^ > 0. Then to satisfy the condition cr^ > 
of the above inequality we define C by 



C = BA-^B'^ ^E, => rank. 



m 



rank E. (13) 



With the above identity the condition that a > is equiv- 
alent to 



BA-^B^ - B^A-^B + E > 0. 



(14) 



We first show that one can choose A > and B and 
^ > such that rank £; = 1 and Eq. (fni) hold. For this 



purpose, we will see that it is enough to find A > and 
B such that G := BA'^B^^ -B'^A'^B has exactly m - 1 
strictly positive eigenvalues and one negative eigenvalue. 

Let F > given. Then F has the spectral decompo- 
sition F = /7A/7^, where U is unitary and A > is a 
diagonal matrix with the diagonal entries equal to the 
nonnegative eigenvalues of F. Choose B = UA^ and 
A = I. Then G = F — A. We claim that we can choose 
F such that G has m — 1 positive eigenvalues and one 
negative eigenvalue. 

Fix H = [hij] G H2m with zero diagonal, i.e ha = 
for i = 1, . . . , n, and a diagonal D = diag((ii, . . . , dm)- 
Assume that di > . . . > dm > 0. Choose t :^ 1 and 
consider H{t) = tD ^ H = t{D ^ \H). Set z = | 
and recall that D{z) = D -\- zH has analytic eigenvalues 
for small z. Since the eigenvalues of D are simple, and 
Dei = diGi^ where e^ = (^li, . . . , ^m)^ it follows that 
the eigenvalues Ai(z), . . . , Xm{z) of D{z) have the Taylor 
expansion 

Xi{z) =di + 0(^^), i = l,...,m 

since H has zero diagonal, (see e.g. p^). By choosing 
F = H{t),t > 1 we deduce that G{t) := H{t) - A{t) = 
H^0{\). 

It is left to show that there exist hermitian H with 
zero diagonal entries and m — 1 positive eigenvalues. Let 
Ai > ... > A^. It is known (Schur's theorem, e.g. J3l 
(5.5.8)]) that the sequence (Ai,...,Am) must majorize 
the sequence of the diagonal entries (0, . . . , 0) of H. 



i=l 



A^, > ^0 = 0,r = l,...,m-l, X]Ai = ^0 = 0. 
1=1 i=i i=i 

(15) 
Furthermore, if Ai > . . . > A^ satisfies the above condi- 
tions, then there exists a real symmetric matrix H with 
zero diagonal and the eigenvalues Ai, . . . , \m (see Theo- 
rem 4.3.32 in [Ml). Choose Ai > . . . > A^_i > and 



\ — _ Y^^--L \ 

H with zero diagonal and m 
values. Hence for t ^ 1 G{t) has m 



Then (15) holds. Thus there exists 

1 strictly positive eigen- 

1 strictly positive 



eigenvalues. Choose to ^ 1 and set G = G{to). Let 
G\u) = Xm\u)^ where Am < 0. Let Eq = —2Xm\u){u\. So 
G^ Eo> and rank Eq = 1. For fc > 1 let ^^i > such 
that rank {Eq -^ Ei) = k. Then E = Eq ^ Ei. D 

Note that from Theorem |3] it follows that we can 
rewrite the expression of the relative entropy of p sim- 
ilar to the formula (7) of Ref. [12 . That is, 

Er{p) = Tr(plogp) - Tr(aloga) +xTr(L-i((/))loga). 

(16) 
From the theorem above it follows that for the case 
m = 2, if <j > then a^ can have at most one zero 
eigenvalue. Hence, for this case (j) is unique as pointed 
out in [12]. For m = 3 it follows from the theorem that 
there exists a > such that a^ has two independent 
eigenstates corresponding to zero eigenvalue. Here is an 
example of such a state a of the form (12) 



zero eigenvector. For the singular a we work below with 
the basis where a is diagonal 



a = diag(5i,...,5n) , 



(17) 



where si > . . . > Sr > = s^+i = . . . = 5^ and 1 < r < 
n. Here r = rank a < n since a is singular. Note that in 
this basis p has the following block diagonal form 



Pii 




, where pu e H^,+,i. 



(18) 



With this eigen-basis of a, we define the matrices 
T(cr), S{a) on the support of a just as in Eq. ([6|, and 
zero outside the support (i.e. the last n — r rows and 
columns of T{a), S{a) are set to zero). Note that with 
this definition 

T{a)oS{a) = S{a) oT{a) = P, = diag(l, . . . , 1, 0, . . . , 0), 



a = 



229 



1 





1 


1 


1 





1 


100 


6 





8 






6 


8 ] 




















46 


60 


60 


80 J 



VI. TENSOR PRODUCTS 

We now show briefiy how to extend the results pre- 
sented in this paper to tensor product of separable states. 
For this purpose we denote by V"^ and V^ the set of sepa- 
rable states in Alice's lab and Bob's lab, respectively. We 
also denote by P^^ the set of separable states of the com- 
posite system. First observe that if aa G dV"^ then for 
any separable state (j[, G V^ the state (7^ (j[, G dV^^ . 
Furthermore, let (j)a G H^ be a supporting hyperplane 
of V^ at (Ja of the form given in Eq. Q and (10). Let 
(j)'^ G H^/, which is nonnegative on P^, i.e. Ti{^^'^) > 
for all a^ G V^ (i.e. ^J, is an entanglement witness in 
Bob's lab). Assume the normalization Tr((0J,)^) = 1. 
Then it is straightforward to show that (j) := (pa^^^b sat- 
isfies Eq. ([2| and ([lo| for any a' G P^^ and a = a^ (8) cr[, . 



Assume first that <Ja > 0, crj, > 0. Then we can use = 
^a^^5 in the formula of Eq. (|3| to find the corresponding 
entangled state p G Hnn',+,i- If cr^ > and a^ is singular, 
we can still use the formula in Eq. (|3|, where (/)5 > on 
V^ and (/)J,x = if a^x = 0. If (7^ is singular then we 
can use the formula given in the next section. 



VII. THE CASE OF SINGULAR CSS 

If the entangled state p is not full rank then the CSS 
a can be singular (i.e. not full rank). More precisely, if 
X is an eigenvector of a corresponding to zero eigenvalue 
then X must also be an eigenvector of p corresponding to 



where Pa the projection to the support of a. Define now 
the linear operators L^r, i^^ : H^ ^ H^ 

L,{0:=T{a)o^, lUO := S{a) o ^. (19) 

Then L^r and L'l are selfadjoint and 

LaLi = LlLa = Pa. (20) 

Note that Lj. is the Moore-Penrose inverse of L^r, and 
that if cr > then Ll = L~^. Note also that L^(cr) = P^. 
With the above definition for L^r, Eq. (p| can be gen- 
eralized to the singular case: 

Lemma 6. Let a G H^^^ be a nonzero singular matrix. 
Let ^ G Hn be positive on the eigenvector subspace of 
a corresponding to the zero eigenvalue, fx'^'^x > if 
(jx = and x 7^ 0.^ Assume that p G H^ is nonzero and 
px = i/ <jx = 0. Then there exists e > such that for 
any t G (0, e) the following hold. 



Tr(plog(a+tO) 



Tr(ploga)+tTr(pL,(0)+O(t2| logtl). 

(21) 



Proof. Without a loss of generality we may assume that <j 
and p of the form ( p!7[ ) and ([l8|. (However we do not need 
the assumption that pu > 0.) Then there exists £ > 
such that for cF{t) > for t G (0,6). Rellich's theorem 
yields that the eigenvalues and the eigenvectors cF{t) := 
(j-ht^ are analytic in t for \t\ < e. Let si(t), . . . ^Sn{t) be 
the analytic eigenvalues of a{t) such that 



Si{t) = Si -^bit 



Em^ 



i = 1, . . . ,n. 



i=2 



The positivity assumption on ^ implies that bi > for 
i = r + 1, . . . ,n. (Note that s^ = for z > r.) Rel- 
lich's theorem also claims that the eigenvectors of o-{t) 
can parameterized analytically. So there exists a unitary 



U{t)^t e [0,£), depending analytically on t, for |t| < £, 
such that the following conditions hold. 

a{t) = U(t) diag(si(t), . . . , Sn{t))U(t)^ 

CO 

U{t)U\t) = In, U{t) = Y.^Wj, Uo = In. 

j=o 
Hence 

log(a + to = U{t) diag(logsi(t), . . ., log Sn{t))U^{t) 
= U(t) diag(0, . . . ,0,log5,+i(t), . . Aog{sn{t)))U\t) 
+ U(t) diag(log5i(t), . . . ,log5,(t),0, . . .,^)U\t). 

Note that the last term in this expression in analytic it t 
for \t\ < e. Clearly, logSi{t) = log{bit) + analytic term. 
Observe next that 



Note also that Ti P^ja' < 1 and therefore for any x G 

[0, 1], '^(x) = xi/j -\- {1 — x)Pa- also satisfies the same con- 
dition Tr('0(x)<j') < Tr('^(x)<j) = 1. We now ready to 
prove the main theorem for the case of singular CSS. 

Theorem 7. Let a G dV he a singular matrix in the 
boundary ofV, and let p G H^,+,i he an entangled state. 
Then, 

Er{p) = S{p\\ct) {i.e. a is the CSS of p) 

if and only if p is of the form 

p{x, a) = (1 - x)a + xLl (V^) , < x^ax < 1. (23) 

Here Xmax < 1 is the maximum value of x not greater 
than 1 such that p{x^a) G H^^+^i. The supporting hyper- 
plane of V at a is represented by ip ^ P^, that satisfies 
the conditions in Eq. pM) and is zero outside the support 
of a (i.e. if CF^ = then iIj^ = 0). 



Uit) diag(0, . . . ,0,logs,+i(t), . . .\og{sn(t)))U\t) 
= diag(0, . . . , 0, log s^+i(t), . . . \og{sn{t))) 
+ tUi diag(0, . . . , 0, log Sr+i(t), . . . log{sn{t))) 
+ t diag(0, . . . , 0, log5,+i(t), . . . \og{sn{t)))Ul + 0{t'\\ogt\). Tr(4(V^)) = Tr(4(V^)P,) = Tr(V^4(P,)) = Tr(V^a) = 1 



Proof. Suppose first that p = p{x,a). We want to prove 
that <j is the CSS of p. First, observe that 



Using the standard fact that Ti XY = TiYX and the 
form of p given by (p^Sl) we deduce that 



TT{pU{t) diag(0, . . . , 0, log s,+i(t), . . . \og{sn{t)))U^{t)) 
= 0{t'\\ogt\). 

Hence 

Tr(plog(a + tO) 

= TT{pU{t) diag(log5i(t), . . . ,log5,(t),0, . . . ,0)/7t(t)) 

^0{t^\\ogt\). 

Similar expansion result hold when we replace a by a 
a diagonal a > as in the beginning of this section. 
Combine these results to deduce the validity of (21 ). D 



Hence p(x, a) G H^^+^i for x G [0,Xmax], where we 
assume that x^ax < 1- It is left to show that 
Tr(p(x,cr)logcr') < Tr(p(x,cr)logcr) for any a' G V. 
From the continuity argument, it is enough to show this 
inequality for all a' > 0. Let a{t) = (1 — t)<j + a't. Let 
f{t) = Tr(p(x, cr) log cr(t)). Using the equality (21), sim- 
ilar to Eq. (11), we get for p = p{x^a): 



From the lemma above it follows that Proposition [T] 
holds true also for singular p and singular CSS a. To 
see that, let a be singular CSS of an entangled state 
p, and suppose first that a^ > 0. Define also a{t) = 
(1 - t)a -^ta' = CF^ t{cr' - a). Note that ^ := a' - a sat- 
isfies the assumptions of Lemma [6J Thus, the arguments 
in Proposition IT] yield that Tt pL(j{a') < TrpLcr(cr) = 1. 
Using the continuity argument, we deduce that this in- 
equality hold for any a^ eV. 

In Eq. ([2]) we defined the supporting hyperplane in 
terms of Hermitian matrix (p satisfying that Tiacj) = 0. 
As we will see below, it will be more convenient to rep- 
resent the supporting hyperplane of P at a in terms of 
i/j = I — (j). That is, the supporting hyperplane will be 
described by the linear functional ^ : H^ -^ R, defined 
by ^(0 = Tr('0^), where '^ satisfies: 



/(O) = TT{pL,{a^ - a)) = Tr [{a^ - a)L,{p)] 
= TT[{{l-x)P^^xiP){a'-a)] = 
= {l-x) Tr(a P^) + X Tr(crV) - 1 < , (24) 

where we have used that Lcr{p) = (1 — x)Pcr-\-xtlj and x G 
(0, 1]. (Note that Tra'P^ < Tra' = 1.) Hence f{0) < 0, 
and since f{t) is concave we have /(O) > /(I), which 
implies that Tr [p(x, a) log a] > Tr [p(x, a) log a']. This 
completes the second direction of the theorem. Moreover, 
note that Er{p) = Tr(plogp) — Tr(plog<j) > 0, since 
p ^ (J. Hence p is entangled. 

Assume now that a is a singular CSS of an entangled 
state p' . Without a loss of generality we may assume 



that a and p' of the form (17) and \lm. Hence, from 



Proposition [T] when applied to the singular case (see the 
discussion above), we get 

Tr(j'L,{p')<Tr(jL,{p') = l 

for ah cr' G V. Denote ^j)' = L^{p'). Note that ^l)' ^ P^ 
since p' ^ a. Hence, with this notations the equa- 
tion above reads Tr?/;V < Tr^/^V = 1. Moreover, by 
definition i\)' is zero outside the support of a. Then 
p' = Ll{L^{p')) = Ll{i)). Hence Eq. ^ holds for 
p = p\ip = ip' and X = 1. Note that if we define tjj by 



Tr(V^crO < Tr(V^cr) = 1 , for ah a' G V. (22) 



^ = t^'^{l-t)P^ 



for some t G (0, 1], then 7/; also satisfies the requirements 
of the theorem. By taking L^ on both sides of the equa- 
tion above, we get p = tp'-\-{l — t)(j. This implies Eq. (23) 
X = t. Note that from the first part of the proof we in- 
deed conclude that a is the CSS to p Hence Eq. (23) 



holds for any x G (0, 1]. This completes the proof of the 
theorem. D 

Recall that Theorem |3] claimed that for p > the cor- 
responding CSS is unique. This is no longer true if p 
semi-positive definite [12 . The reason for that is quite 
simple. 

Theorem 8. Let p G Hn,+,i and assume that p is sin- 
gular. Denote fp{£,) = Trplog^ for ^ G H^^^. Then 
/p(0 > ~^^ f^^ i ^ ^n,-\- 'if tt^^ only if one of the fol- 
lowing condition holds. 

1. ^>0. 

2. pyi = if ^yi = 0. Equivalently, assume that p G 
H^^+ is in the block diagonal form (18) where pn > 
0. There exists a unitary matrix U of order n — 
r such that ^ = diag(/^, U)^ diag(^2, 0) diag(/^, U), 
where < ^2 ^ ^p,-\- ^^^ P ^ [r, n — 1] . 

Denote by Hn,+ (p) the set of all ^ G Hn,+ such that 
fpiO ^ —00. Then Hn,+ (p) is a convex set. The func- 
tion fp : H/e,+ (p) ^ R 25 concave but not strictly concave. 

Proof. Clearly, if ^ > then fp{C) > —00. More pre- 
cisely. Assume that xi , . . . , x^ is an orthonormal system 
of eigenvectors of ^ with the corresponding eigenvalues 
xi > ... > Xn > 0. Then fp{^) = ^^_^(logXi)xJpXi. 
Using the continuity argument we deduce that this for- 
mula remains valid for ^ semipositive definite. Hence 
/p(0 > — 00 if and only if x'^'px = for each eigenvector 
X in the null space of ^. Since p G Hn,+ it follows that 
x'f'px = 4=^ px = 0. This proves the first part of 



|2j The second part of [2] follows straightforward from this 
condition. 

We now show that Hn.+(p) is a convex set. Let ( G 
H^^+ and s > 0. Then g{s) = /p(C + ^^k) is a strictly 
increasing function on (0,oo). (Choose an eigenbase of 
(.) Assume that ^, r? G H^,+ (p). Then for any 5 > 0,t G 
(0, 1), the concavity of log (" on ^ > yields that log(t(^ + 
sh) + (1 - t)(7/ + sh)) > t log(^ + sh) + (1 - t) \og{r] + 
sik), which implies fp{t{£, + sh) + (1 - t){r] + sik)) > 
tfpii + sh) + (1 - t)fp{'n ^slk). Letting s\0 and using 
the assumption that fp{£,),fp{r]) > —00 we deduce that 
fp{t^ + (1 — t)r]) > —00. Hence H^^+(p) is convex. The 
above arguments show also that fp is a concave function 
on Hn.+(p). 

It is left to show that fp is not strictly concave on 
Hn.+ (p). Let ^ = diag(/3,^2),^ = diag(/3, 7^2), where < 
P G H^,+,^2,^2 ^ Hn-r,+ and ^2 7^ m- Clearly, fp{t£, + 
(1 - t)r]) = Trpn log (3 for ah t G [0, 1]. D 

Corollary 9. Let p G Hn,+,i and assume that p is singu- 
lar. Then the set of CSS to p is a compact convex on the 
boundary of V, which may contain more then one point. 



VIII. CONCLUSIONS 

To conclude, given a state a on the boundary of sepa- 
rable or PPT states, we have found a closed formula for 
all entangled states for which a is a CSS. We have also 
shown that if a is full rank, than it is unique. Quite re- 
markably, our formula holds in all dimensions and for any 
number of parties. As an illustrating example, we have 
analyzed the case of qubit-qudit systems and described 
how to apply the formula for this case. 
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Appendix A: The strong concavity of log A 

Definition 1. For an interval int C M let H^(int) be the 
set of all n X n hermitian matrices whose eigenvalues are in 
int. (Here int can be open, closed, half open, half infinite 
or infinite.) Let / : H^(int) ^ H^ be a continuous func- 
tion. / is called monotone^ strict monotone and strong 
monotone if for any C^A G H^(int) the corresponding 
conditions hold respectively: C > A => f{C) > f{A)^ 
C> A^ f{C) > f{A), C>A^ f{C) > f{A). f is 
called concave^ strict concave and strong concave if for 
any C^A G H^(int) the corresponding conditions hold 
respectively: /((I - s)A + sB) > (1 - s)f{A) + sf{B), 
f{{l-s)A^sB) > {l-s)f{A)^sf{B) if rank {A-B) = n 
and 8 e (0, 1), /((I - s)A + sB) > (1 - s)f{A) + sf{B) 
ifA^Bdindse (0,1). 

A well known result is that the functions ft{A) := 
A^^t G (0,1) and log A are strictly concave and strictly 
monotone on J^^((0, oo)). See [13l ^6.6]. In this section 
we show that log A is strongly concave on Hn{{0,oo)). 
This implies that Tr p log a is strictly concave for a fixed 
p > and all a > 0. Hence the CSS a to an entangled 
p > is strictly positive and unique. We need also to 
consider x'''(log 74)x, where x is a nonzero column vector 
in C^, and A is singular and positive. Then it makes 
sense only to consider only those x G U+(A) C C^, 
where U+ {A) is the subspace spanned by eigenvectors of 
A corresponding to positive eigenvalues. For x G \J-^{A) 
we have that x'''(log74)x > — oo. We also agree that for 
each X G C'^\U+(A) x'f'(logA)x = -oo. Then for each 
X G C"^, the function x'''(logA)x is concave on Hn,+ . we 
agree here that 



t(— oo) = — oo = — oo— oo = 



-oo- 



= R— oo for any t > 0. 



Theorem 10. Let A, B e Hn,+ . Then 

R{A,B) := A^B - {A^ B^ -^B^A^) > 0. (Al) 

Furthermore one has the identity 

((1 - s)A^ + sB^)^ = (1 - s)A ^sB-{l- s)sR{A, B)). 

(A2) 
Hence for t G (0, ^] 

(1 - s)A^ + sB^ < ((1 - s)A + sB 

- (1 - s)sR{A, B)Y < ((1 - s)A + sB)\ (A3) 

(1 - 5) log A + 5 log B < log (1 - s)A + sB 

- (1 - s)sR{A, B)] < log((l - s)A + sB). (A4) 



Proof. Le t A, B G H^,+ and assume that R{A, B) is de- 
fined by (|Al}. We claim that R{A, B) > 0. This is a 
straightforward consequence of the Cauchy-Schwarz and 
the arithmetic-geometric inequalities 

Ix'^A^^^xl < ((x'^Ax)(x'^5x))^ < -(x'^Ax + x'^^x), 

Ix'^^^A^xl < ((x'^5x)(x'^Ax))^ < -(x'^^x + x'^Ax). 
Furthermore, R{A^ B) = if and only if A = B. Clearly 
R{A,A) = 0. Suppose that R{A,B) = 0. Then the 
above arguments yield that we must have the equalities 
in the Cauchy-Schwarz inequalities, and equalities in the 
arithmetic-geometric mean for each x. SoA2x = 52x 
for each x. Hence A^ = B^ ^ A = B. 

A straightforward calculation shows the validity of 



(A2). Hence 



f (1-5)^2 +s52J = {1 - s)A ^ sB - {1 - s)sR{A, B)) 
< {l-s)A^sB for s G (0,1). 

Since A 2 is monoto ne, w e deduce from the above inequal- 



ity the inequality (A3) for t = ^. The inequality (A3) 
for t = 2^ follows by induction. Hence ( A4 ) follows from 



(A3) fort: 



9f • 

IS con- 



Assume that t G (0, ^). By assuming that A 
cave and order preserving we deduce from the above in- 
equality (|A3| for t G (0, ^]. D 



Corollary 11. For each p G H^((0,oo)) the function 
Tr(p log cr) is a strict concave function on Hn((0,oo)). 



Proof Let cr^r] G Hn((0, 00)) and assume that cf ^ rj. 
Then R{r], cr) > 0. Hence for s G (0, 1) 

log((l — s)a -\- ST] — {1 — s)sR{a^ r])) < log((l — s)a + sr]) 

=> (1 — 5) log a -\- S log T] < log((l — s)o- + ST]) 

=> Tr(p((l - 5) log cr + 5 log 7/)) < Tr(plog((l - s)a ^ sr])). 



D 



Obviously the above corollary does not hold if p > 
has at least one zero eigenvalue. If cr, 77 > has same 
eigenvalues and eigenvectors which span the range of p, 
i.e. o-U+(p) = 7/U+(p) = U+(p),cr = 7^|U+(p), then 
Tr(p((l — s) log cr + 5 log 77)) is constant for s G [0, 1]. 



